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Abstract 

The effect of an oscillating electric field normal to a metallic surface may be described 
by an effective potential. This induced potential is calculated using semiclassical variants 
of the random phase approximation (RPA). Results are obtained for both ballistic and dif- 
fusive electron motion, and for two and three dimensional systems. The potential induced 
within the surface causes absorption of energy. The results are applied to the absorption 
of radiation by small metal spheres and discs. They improve upon an earlier treatment 
which used the Thomas-Fermi approximation for the effective potential. 



1. Introduction 



When an external electric field is applied to a metal, polarisation charges are induced 
on the surface. At zero frequency, these charges create an additional potential which 
exactly cancels the applied electric field within the interior of the metal. This paper 
discusses the form of the potential at frequencies uj which are small compared to the 
plasma frequency of the metal, u> p , for the case where the external field is perpendicular 
to the surface of the metal. The potential is calculated self-consistently, using two distinct 
simplified versions of the 'random phase approximation' (RPA) approach. The standard 
random phase approximation is discussed in [1]. The formulation of the 'semiclassical' 
variants of the RPA used in this paper was discussed in [2]. 

This effective potential is not directly measurable, but it does have an influence on 
the electromagnetic response of the surface, in particular on the absorption of radiation. 
The results will be used to address a long-standing problem, concerning the theory for 
absorption of radiation by small conducting particles. Standard electromagnetic theory 
predicts that the absorption coefficient a(uj) of a dispersion of small metal particles is 
proportional to uj 2 , and the coefficients are known for simple geometries such as spheres 
[3], or discs with the electric field vector in the plane of the disc [4]. In the standard 
treatment the conductivity of the metal is assumed to be local, with the current density 
at position r proportional to the electric field at r. This assumption is valid if the motion 
of electrons is diffusive. In very small particles the bulk mean free path may be larger 
than the size of the particle, and in this ballistic case the conductivity must be treated 
as a non-local quantity. The frequency dependence of the absorption coefficient has not 
previously received a satisfactory treatment in the ballistic case. The final result of this 
paper will be expressions for the absorption coefficient of the form 

a(oj) ~ K 3 uj 2 , spheres (1.1a) 

a(oj) ~ K 2 uj , discs . 

In the latter case, the electric field is polarised in the plane of the disc. 

Previous papers discussed the absorption of low frequency radiation in particles with 
ballistic electron motion using a Thomas-Fermi approximation: reference [5] treated spher- 
ical particles in three dimensions, and reference [6] discussed discs with the electric field 
polarised in the conducting plane. The surface may be smooth enough that the electron is 
reflected specularly, in which case the motion is integrable due to conservation of angular 
momentum, or it may be a rough surface resulting in ergodic electron motion. The earlier 
works [5] and [6] consider both integrable and ergodic cases. It was shown that in the inte- 
grable case the absorption coefficient is the sum of contributions from resonant absorption 
by electrons with angles of incidence 9 satisfying the condition 

(vy / au)){nn ± 6 n ) = sm9 n (1-2) 

where the integer n labels the resonance, vp is the Fermi velocity, and a is the radius. 
There are no resonances below the frequency uj c = vp/a, so that there is no absorption 
when u < uj c . As the frequency increases, more resonances contribute, and it was shown 
that when uj ^> u c , the frequency-averaged behaviour of the sum of these resonances is of 
the form (1.1). In the case of ergodic electron motion there are no resonances, but (1.1) 
continues to apply (although the coefficients K 2 and K 3 are different). The prediction that 
the absorption coefficient is proportional to frequency in the two dimensional case was a 
surprising result. 



In [5] and [6] it was also shown that when to ^> u c , the frequency-averaged absorption 
can be obtained correctly by treating the collisions of the electron with the surface as if 
they are independent events. The present paper adopts this simplification, which avoids 
summing over contributions from a large number of resonances. 

Reference [2] gave a comprehensive discussion of the equations underlying the treat- 
ment of absorption of radiation, and concluded that the Thomas-Fermi approximation is 
not sufficient, even when the electron motion is ballistic. The form of the potential was 
shown to be very different from the Thomas- Fermi approximation when u ^> u c . The con- 
clusions of references [5] and [6] should therefore be re-evaluated, particularly the prediction 
that a(u>) ~ lo in the two-dimensional case. This paper treats absorption by particles with 
ballistic electron motion using two distinct approximations to the 'random phase approx- 
imation' (RPA), instead of the simpler Thomas- Fermi approximation used in [5] and [6]. 
The conclusions are consistent with (1.1) (but different values for the coefficients K 3 and 
K 2 are obtained). 

Two simplifications of the RPA approach are discussed in section 2. The first will be 
termed the 'image source approximation', and leads to a slight simplification of the RPA 
equations. It is valid in the neighbourhood of a flat surface, when the frequency is small 
compared to the plasma frequency. The second approach, previously discussed in [2], will 
be termed the 'semiclassicaP RPA method. It makes more radical assumptions, and leads 
to simpler equations. The results of both approaches are approximate, rather than being 
the leading term of an asymptotic theory. Section 3 obtains the self-consistent effective 
potential which describes the electric field within the surface. First this is calculated 
using the semiclassical RPA method for both two and three dimensional systems, and then 
using the image source approximation for the simpler case of a three dimensional system. 
Section 4 considers the transfer of energy to an electron colliding with the surface. Again, 
the calculation is carried out first for the semiclassical RPA method, and then using the 
image source approximation in the three dimensional case. The rate of absorption of energy 
is calculated in section 5, for both spheres and discs. Section 6 is a brief discussion of the 
validity of the results. 

This work uses a free electron model for the conduction electrons, as described in 
[5]. The symbol e will denote the magnitude of the electron charge, and the potential 
<fi will denote the potential energy of an electron (rather than the electostatic potential). 
Following common practice equations will be written as equalities, despite the fact that 
most of them are approximate relations. References [2] and [5] contain ample discussion 
of and references to recent literature. 

2. Two semiclassical RPA methods 

2.1 The RPA equations 

A polarisable medium is perturbed by applying a time-dependent external potential, 
which is specified in the frequency domain by a function </> e xt(r, lo). The motion of electrons 
within the medium may be approximated by an independent particle effective Hamilto- 
nian, which contains an effective potential 0(r, to). By analogy with Dirac notation, these 
potentials can be notated as function space vectors |0 ex t) an d \<j>) respectively. The effec- 
tive potential must take account of the fact that the medium is polarised by the externally 
applied field, so that the effective potential is the sum of the external potential and the 
potential |0 po i) generated by the polarisation charge \p) through the action of the Coulomb 



operator U: 

pol (r, u ) = (r|0 pol ) = (r\U\p) = ^- q J dv'y^p(v\ u) . (2.1) 

The induced charge density \p) is obtained from the effective potential |0) by multiplication 
by the polarisability operator II(cu), i.e. \p) = Yl(uj)\(p). Writing |0) = |</> ex t) + \4>poi) gives 
a single equation which should be solved for the effective potential: the RPA equation is 

|0) = |0ext) + ^n(o;)|0) . (2.2) 

A thorough discussion of the RPA method is given in [1] . 
2.2 The image source approximation 

The polarisability TI(uj) is related to the spatial probability propagator by 

tl(uj) = ev[i + iujP(uj)] (2.3) 

where (r|P(o;)|r / ) is the Fourier transform of the probability P(r,r',t) of an electron ini- 
tially at r' being located at r after time t. A semiclassical derivation of this relation is 
given in [2]. 

In the vicinity of a flat surface, the coordinate space representation of the polarisability 
may be approximated as follows: 

U(r,r',oj) = (r|n(w)|r') = ev[8{r - r') + iuP{r - r',oj) + iojP(r - r' R ,oj)] . (2.4) 

where P(r — r',oj) is the Fourier transform of the free space propagator P(r — r',t), and 
tr is the reflction of the point r in the plane of the boundary. 

Specialising to the case where the boundary is a flat plane at z = 0, and where the the 
potential (f>(z) depends only on the distance from the boundary, the equation \p) = Il(c<;)|(/>) 
may be written 

p(z) = eu(j)(z) + ieuu J dx' J dy' J dz [P((x' \y' , z — z'),u) + P((x\y' ', z + z'),u)~\(j)(z') . 

(2.5) 

The integrals will be taken to be over all space, with 4>{z) = for z < 0. Equation (2.5) 
will be termed the 'image source approximation'. 

In view of the isotropy of the free space propagator, its Fourier transform P(k, u) is 
a function of the magnitude k = |k| of the wavevector: 

p(k,u) = P(k,u>) = J dr expfik.r] P(r,u) , k = |k| . (2.6) 

The function p(k, uj) can be determined by Fourier transformation of a semiclassical 
approximation of the position representation of the free space propagator, valid when 
k /cf- In [2] it was shown that the resulting Fourier representation of the free space 
polarisability may be expressed in terms of a single scaling variable A: 

IL(k,u)) =ei/[l + iwp(k,u))] =eu[l+g(X)] . (2.7) 



The dimensionless variable A takes different forms for ballistic or diffusive electron dynam 
ics 

kv-F/u , ballistic 



A = ^ li^' (2.8) 
[ky/D/oj , diffusive 

where D is the diffusion constant. In the ballistic case the form of the function g(k) 
depends upon the dimensionality of space. In two dimensions: 

^ 2(A) = \ i(A2-l)-i/ 2 , A>1 (2 - 9) 



and in three dimensions 

A + l 



f3(A) = -7^ lo Se 



A- 1 



+ He (A -l) (2.10) 



(where 0(x) is the Heaviside function, with increasing unit step at x = 0). The limiting 
forms of fir (A) for small and large argument are informative: 

S<J (A) = -(1 + iA 2 ) + 0(A 3 ) (2.11) 

In both two and three dimensional ballistic motion, <7d(A) has an integrable divergence at 
A = 1. In the case of diffusive electron motion, the form of the function g(X) is the same 
in two and three dimensions: 

0d(A) = • (2.13) 

2.3 The semiclassical RPA approximation 

Formally, the effective potential may be determined simply by inverting (2.2): |0) = 
[I — UIl(uj)]~ 1 \(/) ex t) ■ This approach can be implemented numerically by expanding func- 
tions in a suitable basis set and inverting matrices numerically. It is desirable to have a 
simpler approach which allows further analytical progress. In [2], a 'semiclassical' variant 
of the RPA approach was introduced. It is assumed that (provided w < w p ) the polar- 
isation charge \p) is that which would be predicted by classical electrostatic theory (the 
justification for this is discussed in [2]). Under this assumption, the effective potential 
satisfies a much simpler equation 

|p cl )=n(u;)|0) (2.14) 

where \p c \) is the classical charge distribution function resulting from a static external 
field, which is assumed to be known. The task of determining the effective potential is 
then reduced to the simpler task of inverting ft (a;). Calculation of |p c i) is still a difficult 
problem, but solutions are obtained in various geometries in textbooks such as [3] . 

The form of the surface charge density is different in two and three dimensions. In 
three dimensions the charge density is confined to a narrow layer at the surface of the 
conductor, with thickness equal to the Thomas-Fermi screening length. This can be ap- 
proximated by a delta-function distribution, a tiny distance e inside the surface: 



p(z) =qS(z-e) 



(2.15) 



where q is the surface charge density induced by the externally applied fields. In (2.15) 
the coefficient q depends upon the position s on the surface of the conducting particle. 
In the case of a two-dimensional conductor in three dimesional space, there is an inverse 
square-root divergence of the charge density at the surface of the particle: 



C_ 

7^ 



(2.16) 



The reasons for this behaviour are discussed in [6]. 

3. Calculation of the surface potential 

3.1 Method for solving the semiclassical RPA equation 

Equation (2.5) will now be expressed in a purely one-dimensional form. To this end, 
define F(z, u) as the inverse Fourier transform of p(k, u) 

/oo 
dz exp[ikz] F(z, to) . (3-1) 
-oo 

and note that 

/OO /"OO /"OO 

dx' dy' dz' P((x',y',z- z');u) <j){z') 
-oo J — oo J — oo 

= ~(27)3 / dr "^ z - z ") J dK exp[-iK.r"]P(K;w) 
= -— [ dz" <p(z - z") ( dK P(K, u) 5{K X ) 8{K y ) exp[-iK.r"] 

27T J-oo J 

-I poo poo 

= / dz"cj){z-z") / dk exv[-ikz"]p{k,u>) 

^ J — oo J — oo 

/oo 
dz F(z- z',u)(j){z') . (3.2) 
-oo 

Using this result, the integrals over x' and y' in equation (2.5) may be eliminated, giving 



p(z) = ev 



/oo 
dz [F{z-z ,u) +F(z + z',uj)] (p(z') 
-oo 



(3.3) 



Comparison of (3.3) with (3.1) shows that F(z, u>) is related to the inverse Fourier transform 
of the function #(A), introduced in (2.7). In the ballistic case: 



1 f°° 1 

F(z,u) = / dk exp[— ikz] g{kv-p/u) = — G(zco/vf) 

2-kiuj J x ' W|- 



(3.4) 



where G is the inverse Fourier transform of g. Equation (3.3) can now be written in the 
scaled form 



p(z) = ev 



1 



z — z 



■ z + z 



A 



A 



(3.5) 



where the scale length A is 

A _ f v F /u> , ballistic g ^ 
\ y^D/u , diffusive 

Defining a scaled distance by j = z/A, and a scaled charge density f(x) and potential 

^(x) by 

p(z/A) = evf(Ax) , 0(z) = ^(Ax) (3.7) 
equation (3.5) may be expressed in the dimensionless form 

POO 

f(x) = ip(x)+ dx' [G(x-x') + G(x + x')]iP(x') . (3.8) 
Jo 

Now consider how to solve equation (3.8) for the scaled potential ij)(x), given the scaled 
charge density f(x). The function f(x) is defined only for x > 0. Also, the behaviour of 
tp(x) for x < is irrelevant to the form of f(x) in the region where this is defined. We 
may define a symmetric extension of the function f(x): 

fs(x) = f(\x\) ■ (3.9) 
Consider a function ip a (x) which satisfies 

/oo 
dx' G{x-x')il) a {x') . (3.10) 
-oo 

This function must be symmetric: ip a (x) = ip s (—x). This symmetric solution satisfies 
the same equation as ip(x) (equation (3.8)) for x < 0. We therefore solve the simpler 
equation (3.10), and drop the subscript s labelling the solution. The Fourier transform of 
the solution is 

^ = irSJj ■ (3 - n) 

The solution may also be expressed as a convolution 

/oo 

dx' K{x - x') f 8 {x') (3.12) 
-oo 

where K{x) is the inverse Fourier transform of (1 + g(k))~ 1 . 

These expressions require the form of the function f s (x). Referring to (2.15) (and 
taking e — > 0), in three dimensions the Fourier tansform of the scaled and symmetrised 
charge density is, in the ballistic case 

/.(A:) = . (3.13) 

The factor of 2 appears in (3.13) because both the delta function and its symmetric image 
contribute. From (2.16), in two dimensions the analogous quantity is 

UK) = ^pZ (3 . 14 ) 



3.2 Solution of the image source approximation 



It is instructive to compare the solution of the 'semiclassicaP RPA equation, (2.14), 
with that of the image source approximation, (2.5). This is difficult in the two-dimensional 
case, but quite straightforward in the case of a flat metallic surface in three dimensions. 
The charge density is still assumed to be a function of z alone, and is related to the effective 
potential <f>(z) by Poisson's equation: 

i£| = *« (3.15) 

e dz z eo 

so that the charge density in (3.5) may be replaced by a term proportional to the second 
derivative of <j>(z). (This simplification is not possible in two dimensions, because it that 
case p also depends upon the second derivative with respect to the coordinate perpendic- 
ular to the conducting plane.) The resulting equation can be expressed in a scaled form 
analogous to (3.8), and is then transformed into an equation for the symmetrised potential, 
corresponding to (3.10). Here it is necessary to note that the symmetrised potential can 
have a discontinuity at z = 0, without the charge density having a singularity there. The 
electric field approaches a constant value as z — > oo, and provided w <C w p this field is 
much smaller than the externally applied field. If the internal field is neglected, the po- 
tential <f>(z) may be assumed to approach zero as z — > oo, and as z — > + the slope d(f>/dz 
approaches —eq/eo, where q is the integral of p(z), i.e. the total charge per unit area bound 
to the surface. The symmetrised potential therefore has a discontinuity of slope equal to 
—2qe/eo at z = 0. The scaled potential ip 3 (x) satisfies 



e 2 z/A 2 dx 2 



= ip s (x)+ f dx' G(x-x')^Jx')- ^-8(x) . (3.16) 
J-oc euA 



This version of the RPA equation can also be solved directly by a Fourier transform ap- 
proach. Noting that for ballistic electron motion in three dimensions eo/e 2 z/A 2 = cu 2 /3o; 2 , 
the Fourier transform of the image source approximation solution is 

^{k) = ^L * (3.17) 

This solution is indeed very close to the semiclassical RPA solution when uj <C uj p . 

It should be emphasised that these forms for the effective potential are applicable only 
when the scale length A is small compared to the characteristic dimension a of the system. 
In the low frequency limit, uj <C uj c the scale size A exceeds the size of the particle. In this 
limit a Thomas-Fermi approximation should be used for the effective potential. In this 
case, the form for the Fourier transform of the effective potential is obtained from (3.11) 
or (3.17) by replacing the function g{k) by zero. 

3.3 Components of the effective surface potential 

Reference [4], which considered the case where electron motion is diffusive, suggested 
writing the effective potential as the sum of two terms, namely a 'static' potential, which 
is given by the Thomas-Fermi approximation, and a 'dynamic' potential, which is propor- 
tional to frequency, and which is required to move the static polarisation into place when 
the external electric field changes. Reference [2] showed that an additional component 
must be present in the ballistic case. It is interesting to see how these components are 



represented in equations (3.11), (3.13) and (3.17). The discussion will be restricted to the 
three dimensional case. 

Thomas-Fermi theory predicts that the potential is proportional to the charge density: 
(j)(z) = p[z)jev. This is equivalent to setting g(k) = in (3.11) or (3.17). The resulting 
potential will be termed the 'static' potential, <f> a ta,t(z). I n the case of the semiclassical 
RPA method, it is simply a delta function localised at the surface, and in the case of the 
image source approximation it decays rapidly as a function of distance from the boundary, 
with decay length A s = y^eo/e 2 iy. 

In the vicinity of a surface which accumulates a polarisation charge q, the internal 
electric field has magnitude S[ nt = iuq/cr(u) where cr(u) is the bulk conductivity. Reference 
[2] showed that this expression for the internal field is also applicable in the ballistic case, 
when the distance from the boundary is greater than A = v-p/oj. In the ballistic case, the 
bulk conductivity is c(oj) = Ne 2 /imui, where N is the density of conduction electrons, and 
in the diffusive case a = e 2 uD for frequencies small compared to the collision rate. There 
is therefore a component of the symmetrised potential which is proportional to \x\. In the 
ballistic case, this component of the scaled potential is 



V'dyn(^) 



ieooqA , 



x 



3quJ 



2evvY 



\x\ 



(3.18) 



The generalised Fourier transform of \x\ is —2/k 2 . Consistency with (3.18) therefore re- 
quires that the Fourier transform of ip s (x) approaches 3qu / 'ew-pk as k — > 0. Using 
equation (2.11) to evaluate the limit of (3.11) or (3.17) as k — > verifies this relation. 

The full effective potential may be written as the sum of </> s tat (z) , (f>dyn(z), and an 
additional term 'surface' term <j> 8VLT f(z), which decays with a characteristic length scale 
A = v-p/u. The three components of the effective potential are illustrated schematically in 
figure 1. The surface potential for the three dimensional ballistic case is 



VWf(aO 



2qu 
euvp J 



dk exp(ikx) 



j_ 

2k 



log e 



fc+1 
fc-1 



+ £6(fc -!) + §(£): 



~ l+ k 2 



(3-19) 

This expression diverges logarithmically as x — > when — = 0. 

UJp 

In the case of diffusive electron motion it was shown in [7] that the effective potential 
is the sum of the static and dynamic components only. It is instructive to see how this 
conclusion is confirmed using the results of the present paper. In the diffusive case, it 
follows from (2.15) that the Fourier transform of the charge density is p s (k) = 2q/eu. 
Using (2.13) and (3.11), the Fourier transform of the effective potential is therefore 



k{k) 



2q 1 _2q 

evl + g D (y/D/u)k) eu 



iuj 

m 2 



(3.20) 



The first term in the final bracket Fourier transforms into a delta function, and therefore 
represents the static potential. The term proportional to u/k 2 Fourier transforms to a 
term proportional to uj\z\, and represents the dynamic potential. The additional surface 
potential component is therefore absent in the diffusive case. 



4. Energy transferred on collision with surface 



The objective is to calculate the energy transferred to an electron from the externally 
applied electromagnetic field as it collides with the surface. From this point onwards, 
the discussion is specific to ballistic electron dynamics. The problem will first be treated 
classically, then this will be compared with the results of a quantum mechanical calculation. 

4.1 Classical treatment 

Classically, the energy transfer is determined by separate contributions from collisions 
with the surface of the particle, provided the surface potential is localised at the surface 
of the particle. If the electron collides with the surface at time to, and is in the vicinity of 
the surface for a time At, the energy transferred at this collision is 



6E= dt^(r(t),t) 



-t +At 
't -At 

where r(t) is the trajectory of the electron, and where the potential 0(r, t) is 



(4.1) 



4>(z, t) = Re 



exp(iu;£) ip(z(t)/A) 



(4.2) 



For an electron incident at an angle 9 from the normal to the surface, the distance from 
the surface at time t is 

z{t) = v F cos6\t-t \ . (4.3) 

Using the fact that ip(x) = <f>(xA) is symmetric about x = 0, the energy tranferred is 
therefore given by 



SE = —u Im 



/oo 
dt exp(iotf) ^(^f cos 0(t — t )/A) 
-oo 



(4.4) 



This may be expressed in terms of the Fourier transform of ip: 

-A 



SE = 



V F cos 9 



-Im 



exp(io;to) / dx exp(ixuA/vF cos 9) ip(x) 



-1 

cos 9 



-Im 



exp(io»to)V'(l/ cos6>) 



(4.5) 



Using (3.11), this result may be expressed in terms of the Fourier transform of ^(z) as 
follows: 



SE = 



1 



cos 9 



-Im 



exp(iu>to)- 



/ B (l/cos0) 



(4.6) 



1 + 2 (1/ COS 0). 

Now consider the form of 8E for collision of a ballistic electron with the surface. In 
the two dimensional case, combining (4.6), (2.9) and (3.14) gives 



ev v 



exp(iwto) 



Vcos0(l + icotfl) 



(4.7) 



which may be written 



^ C 2tvu sin0 ^ 
SE = — a ■ Re 



V v F v / cos0 



exp[i(wt + X2(9))] 



(4.8) 



where Xi{^) * s a real- valued phase. In the three dimensional case the analogous expression, 
obtained using (2.10) and (3.12), is 



5E = 



2qu 

euv-p 



Im 



exp(L;to) 



cos#(l + g 3 (l/ cos 0)) 



2qu 
euv-p 



S(6)Re 



exp[i(ujt + x 3 (0))] 



(4.9) 



where 



S(6) 



cos 9 



1-| cos(9[log e (l + cos/9) - log e (l - cos/9)] + § icos6> 



(4.10) 



and Xz{Q) is a phase. Both the two-dimensional and the three-dimensional expressions 
diverge at grazing incidence, where 9 — > tt/2. Physically, this can be interpreted in terms 
of the wavepacket spending a long time in contact with the barrier for a reflection at 
grazing incidence. A quantum mechanical treatment will take account of the fact that the 
wavefunction approaches zero at the surface. This removes the divergence. 

4.2 Quantum treatment 

Now consider how the expression (4.6) must be modified to take account of quantum 
mechanics. The approach will be to consider the collision of a wavepacket with the surface. 
The expectation value of the energy transferred will be calculated: this is 



(5E} = J dt J dz P{z,t)-£(z,t) 



(4.11) 



where P(z,t) = t)\ 2 is the probability density for the electron to be at a distance 
z from the surface. If the wavepacket can be chosen to be sufficiently well localised, the 
energy transfer can be assumed to be equal to this expectation value. On the other hand, if 
the length scale L over which the wavepacket is localised is too large, then the expectation 
value will represent an average over the temporal variation of the potential, and (SE) will 
under-estimate the magnitude of the energy transferred. The criterion is L <C v z /ui, where 
v z = v-p cos 9 is the speed at which a Fermi-surface electron with angle of incidence 9 
approaches the surface. If the energy of the electron is to be well defined (and close to the 
Fermi level), then kpL 1, where kp is the Fermi wavevector. These two inequalities for 
L are compatible provided huj <^ Ep, which is assumed throughout. 

The wavefunction of the wavepacket which collides with the surface at time t = 
may be written in the approximate form 

(z, t) = exp(-i£ F t/ft) [exp(ip z z/h)f(z - v z t) - exp(-ip z z/h)f(z + v z t)] (4.12) 

where f(x) is a symmetric function, which decays rapidly when \x\ ^> L. Dispersion of the 
wavepacket is unimportant, and is ignored in writing (4.12). The function / is normalised, 
and its autoconvolution F is required: 



/oo 
dx' f(x-x')f(x') , F(0) = 1. 
-00 



(4.13) 



Assuming that SE = (5E), and substituting (4.12) into (4.11) gives 

dz / dt -£(z, t) [f 2 (z-v z t)+f(z+v z t)-2 cos(2p z z/h)f(z-v z t)f(z+v z t 

(4.14) 



The integral over time will be performed first. It has been assumed that the envelope 
function varies much more rapidly than the variation in time of the potential (j)(z,t), 
although the potential may vary rapidly as a function of z. Using (4.13), the energy 
transfer may therefore be approximated as follows: 

= J- | dz ^(z, z/v z ) + ^(z, -z/v z ) - 2 cos(2p z z/h)F(2z)^(z, 0) . (4.15) 

If the potential varies sinusoidally in time, such that (f)(z, t) = Re[ex.p(iu>t)(j)(z)] : then (4.15) 
becomes 



5E 



—LO 



V p COS 9 



-Im 



dz exp(iuz / cos 9)<p(z) — F(2z)4>(z) cos(2mv-FZ cos 9 /Ti) 



(4.16) 



The first term of (4.16) is a Fourier transform of 4>(z). Comparison with (3.17) shows that 
this takes the form 

1 



2q 

ev 



l + gs(kv F /u;) + ^ 



(4.17) 



The second term of (4.16) is the Fourier transform of the product F(2z)<p(z) evaluated 
at k = 2mv-F cos 9/fi. This may be obtained by convolution of (4.17) with the Fourier 
transform of F{2z). The function (4.17) has support y 1 e 2 v /eo = 1/A S , where A s is the 
Thomas-Fermi screening length, and has a structure close to k = associated with the 
function g, with a narrower support, uj/v-p = 1/A. The assumptions concerning the support 
L of F(z) imply that A ^> L ^> A s . When calculating the Fourier transform of <p{z)F{2z) 
using the convolution theorem, the structure in (j)(k) associated with the function g is 
suppresed, because the support of F(k) is broader. The function <p(k) remains unchanged 
in other respects, because the support of F(k) is narrower than the overall support of (p(k). 
The quantum mechanical expression for the energy transferred, calculated using the image 
source approximation for the effective potential, can therefore be approximated as follows 



5E = 



-2quJ 



evvp cos 9 



-Im 



1 



1 



1 + <7 3 (/™fA>) i + if£^i C os 2 # 

ve 2 U z 



(4.18) 



This simplifies to 



5E = 



2quJ 
evv-p 



Im 



cos 9 — r<73 ( 1 / cos 9) I cos 9 



exp(iu;to) 



. (r + cos 2 9)(1 + g 3 (l/ cos 9)) 
where the time to of collision with the surface has been inserted, and where 



(4.19) 



4eom 2 Vp 

The constant V was introduced in [5]. It may be expressed in the alternative forms 

r _ 2A F _ 2 5 / 3 r s 
7rao 3 2 / 3 7r 4 / 3 oo 



(4.20) 



(4.21) 



where ao is the effective Bohr radius, and r s is the radius of a sphere containing a single 
electron. Equation (4.19) will also be written in the form 



5E 



2qu 
evvp 



S(9) cos[ut Q + X3 



(4.22) 



where 

cos/9 - r# 3 (l/cos0)/cos0 



(r + cos 2 9)(l+g 3 (l/cos9)) 



(4.23) 



The formulae above should reduce to the Thomas-Fermi theory when the function g is 
replaced by zero, and (4.20) does indeed reduce to equation (3.25) of reference [5] upon 
setting (7 = 0. 

The divergence of (4.10) at glancing incidence is absent in this more sophisticated 
quantum mechanical treatment, because the wavefunction of the electron vanishes at the 
surface, and the electron is therefore unable to be influenced by the potential there. Equa- 
tion (4.19) approaches the prediction from (4.9) and (4.10) in the limit r — > for all values 
of except tt/2, because as V — * the Fermi wavelength becomes small compared to the 
Thomas-Fermi screening length, and the electron is able to penetrate closer to the surface. 

In figure 2 the energy transferred to a reflected electron is plotted as a function of the 
angle of incidence 9, for two different values of T, and for the semiclassical RPA, which 
corresponds to the limit V — > 0. 

5. Rate of absorption of energy 

5.1 General considerations 

The absorption coefficient of a suspension of small particles is determined by the rate 
at which an individual particle absorbs energy. A semiclassical approach developed in 
earlier papers [5,6,2] shows how the rate of absorption of energy by the electron gas may 
be expressed in terms of the variance of the change of energy of a single electron. For 
non-interacting fermions, the rate of change of the total energy of the electron gas may 
be written 

^ = 2VuD E (5.1) 

where V is the volume of the particle, and D E is a diffusion constant for single electron 
energies. The factor of 2 in (5.1) accounts for spin degeneracy, and v is interpreted as the 
density of states per unit volume per spin. The diffusion constant is defined by writing 

(AE 2 (t)} = 2D E t (5.2) 

where AE(t) is the energy transferred to an electron after time t. 

In the context of this paper AE(t) is the sum of the energy 5Ej transferred on collisions 
of the electron with the surface: 



AE(t) = J2*Ej (5-3) 



where the sum runs over all collisions between times and t. The variance in (5.2) is 
defined in terms of a phase space average for electrons at the Fermi energy: 

/a p2/.\\ [da AE\a)d(H(a)-E F ) 

{AE {t)) = Jda S(H(a)-E F ) (5 ' 4) 



where a = (q, p) are phase space coordinates of an electron, and AE(a) is the energy 
transferred to an electron which is initially at a. 



In the ergodic case, 



(AE 2 (t)) = R(5E 2 )t 



(5.5) 



where R is the rate of collision of particles with the boundary. The total distance travelled 
by the electron in time t is vpt = N(L), where (L) is the mean distance travelled between 
each of the N collisions. The rate of collisions is therefore 

R=l*. (5.6) 

5.2 Absorption by conducting discs 

The simpler case in which to evaluate the phase space average (5.4) is for two- 
dimensional discs, and this will be discussed in some detail to illustrate the approach. The 
calculation must be confined to the simpler semiclassical RPA approximation, because the 
more precise image source approximation has not been calculated in the two-dimensional 
case. The integrable case will be considered first, followed by the ergodic case (which 
applies when the surface of the disc is rough). 

For the case of integrable electron motion in a disc, the angular momentum J is a 
conserved quantity, and the coordinates (E, J, to, <A)) are a canonical set, where E is the 
energy, to the time since the most recent collision, and <po is the polar angle of the most 
recent collision. The angular momentum and the period r between collisions are both 
related to the angle of incidence, 9: 

J = mavpsm9, r = — cos6> (5.7) 

where a is the radius of the disc. Also, the charge density induced on a disc by an electric 
field £ in the plane of the disc is 

p(r , v) = i J°^ (5.8) 
7rV a z — r z 

where (p measured from the axis of the electric field. The coefficient C occuring in (2.16) 
is therefore given by 

C(s) = eoSy/a cos ip (5-9) 

7T 

where s = a<p is the distance around the perimeter. In two dimensions, the density of 
states per spin is 

(5 ' 10) 

From (4.8), the total energy transferred to a single electron is 

. . 2nfi 2 l2nuj sin6* ^, s r ,„ Nl 

AE(t) = ^SEj = -—^- rn ==Y. C ^)^ t 3+X2{9)] (5.11) 

me V v f Vcos 9 

J 3 

where tj = jr + t are the times of the collisions with the surface, and the sum runs over 
N ~ t/r values of the index j. The sum in (5.11) is dominated by resonances satisfying 
(1.2). However, in [5] it was shown that when u> ^> u c (5.11) can be approximated by 
assuming that the bounces are independent events. This correctly describes the average 



behaviour of the absorption, but not that of the resonances. The mean squared energy 
transfer at fixed angle of incidence 9 is then estimated to be 



,. ^ T 64nh 4 e 2 S 2 au) sin 2 9 , 9 , . . 9/ ..... 

(AE\t))\ e = N — ^- (cos 2 j )(cos 2 (^ + X2(^))) 

_ 8nh 4 el£ 2 u sin 2 9 _ sin 2 9 

— 5—5 ir-^t = A — t (5.12) 

m z e 2 cos 2 9 cos 2 9 

where the final equality defines A. Now consider how to calculate the phase space average 
in (5.4). The average over ipo has already been performed. It remains to average over J 
and To- The quantity being averaged is independent of To, so that integration over To gives 
a contribution r. The required average is therefore 

2x _fdJr(AE 2 )\ e 



(AE } = 



£ /2 d9 cos 2 9(AE 2 )\ e _ ^ ' dB sin 2 9 
f " /2 d9 cos 2 9 ' f* /2 d9 cos 2 9 ' 



fdJ> 

= M *l /9 Z \ = M ■ ( 5 - 13 ) 



In the integrable case, the rate of absoprtion of energy is dE T /dt = na 2 vA. The final 
expression for the rate of absorption of energy in the integrable case is then 

dE Angela 2 £ 2 u 

— = . (5.14) 

dt me 2 K 1 

It is interesting to note that this expression is independent of the Fermi energy. It may be 
written in the form 

dE T = (eaS) 2 (huj) 
dt ~ 4ithE K [ ' ' 

where Er = me 4 / '16n 2 'e^h 2 is the Rydberg energy, (aeS) is a measure of the energy asso- 
ciated with displacement of an electron across the particle by the electric field, and hu is 
the photon energy. This result differs by a numerical factor from that obtained in [6] using 
the Thomas-Fermi approximation for the effective potential. 

If the edge of the disc is rough so that the motion is ergodic, the rate of absorption of 
energy is calculated from (5.5) and (5.6). Using (4.8), (5.9) and (5.10), the mean squared 
value of the energy transferred at a single collision is 

m I e z v-p \ cos 9 / 
The average over the angle of incidence is 



= J dJrsm 2 9/cos9 j j dJr = d9 cos 9 sin 2 9 j 1 d9 cos 2 9 = ^- . 

(lit) 



Also, 

r/2 / /•tt/2 



// f f n/ I f n/ 16a 

dJr 2 / JdJT = 2aJ d9 cos 3 9 / J d9 cos 2 9 = — . (5.18) 



In the ergodic case, the rate of absorption is 



dE T 2nh 2 ela 2 E 2 uj 

— - = . 5.19 

dt me 2 K 1 

This differs by a factor of | from the integrable case, given by (5.14). 

5.3 Absorption by conducting spheres 

The calculation of the energy absorbed proceeds by analogy with that for conducting 
discs. The first step is to specify a convenient set of phase-space coordinates. In the 
case of ballistic and specularly reflected electrons moving in a spherical enclosure, angular 
momentum is a conserved quantity. The following variables can be used to specify the 
phase-space coordinates of the electron: its energy E, angular momentum vector J, and 
two angle variables, <p and <p'. In [5], it was shown that the measure da for canonical 
coordinates is given by 

da = T - dE dip dip' dJ . (5.20) 

rJ 

The energy transferred to an electron bouncing at the surface is given by equations (4.19) 
or (4.22). The charge density at the surface is that given by classical electrostatics, namely 

q = 3eo£cosx (5.21) 

where S is the amplitude of the external electric field, and x is the polar angle of the point 
on the surface measured from the direction of the electric field. 

The rate of absorption of energy is calculated using (5.1) and (5.2). In the case of 
integrable motion (specular reflection at the surface) , the mean squared energy transferred 
to the electron (AE 2 ) is obtained by first averaging 5E 2 at fixed angular momentum J, 
and then integrating with respect to angular momentum. The variance of SE at fixed J is 

a 2 2c2 

{8E 2 )\j = ^lf^ CO s\ S 2 {9) = Acos 2 Xo S 2 (0) (5.22) 

where %o is the angle between the angular momentum vector J and the direction of the 
external electric field, and the final equality defines A. The number of collisions per unit 
time is 1/r. Treating the collisions as if they are independent events gives 

(A 2 E(t)) = At J dJ ~cos 2 xoS 2 (9) j J dJj 



At I dJ I dxo2ivJ 2 sm X o jS 2 (6) cos 2 xo j J dJ4ivJ 2 j 



o 2 2c2 /-tt/2 / r 7r/2 

' j d9 cos 6 sin eS 2 {9) / j d0 cos 2 6>sin6> (5.23) 



2e 2 v 2 avY Jo 

The rate of absorption of energy by the spherical particle is therefore 



dt e 2 WY 

where 



dE _ ««8«W (5.24) 



f°° 1 \l-Tx 2 g(x)\ 2 
F(T) = / dx 1 ^ 7 . 5.25 

h ^(i + rx 2 ) 2 |i + (7(x)| 



The integral ^(T) diverges logarithmically as V — > 0: for small r, T ~ K — |log e (r), 
where X is a constant. In the limit r — > oo, JF approaches a finite limit. Some values of 
^-"(r) obtained by numerical integration are given in table 1. A clearer understanding of 
(5.24) is obtained by expressing it in terms of ratios of energies: two equivalent forms are 

— E = &Lfr*g)'(MV (r) = ^.MW^ (5 . 26) 

In the case where the surface of the spherical particle is rough, the electron motion 
is ergodic, and the rate of absorption is calculated via the microcanonical average of 5E 2 , 
using (5.5) and (5.6). A complicating feature is that the charge density concentrates 
on prominences of a rough surface. In [5] a simplified model was discussed, in which a 
fraction rj of the surface comprises high plateaus, with the the charge density is increased 
by a factor I/77, and the remainder of the surface is un-charged. According to this model, 
5E 2 is increased by a factor of l/rf for a fraction 77 of collisions. The required average is 

(5E 2 ) = ^ JdJ 1 cos 2 X S\e) / JdJj 

= ^ J dJ j d X 2ii J 2 sin X j cos 2 X S 2 (0) j J dJAnJ 2 ^ 
2A r 17 / 2 

= — d9 cos 2 6 sin 9S 2 {9) (5.27) 
V Jo 

where A is the factor defined in (5.22). The average distance between bounces is 

(L) = v F (r) = Y- ( 5 - 28 ) 
The rate of absorption is then found to be 

^ = ig^w 

dt r\e l vv? 

where 

g(r)= r rf 4 11 -w (5 . 30) 

A ^ (i + rx 2 ) 2 |(i + (7(x)| 

The function (?(r) appraoches finite limits as T — > ±00. Values obtained by numerical 
integration are given in table 1. 



6. Concluding remarks 

The principal new results in this paper were the solutions of the simplified RPA 
equations for the effective potential <p(z) (section 3), and the calculation of the energy 
transferred to an electron rebounding from the surface (section 4). They improve upon 
the earlier analysis in [5] and [6], which used a Thomas- Fermi approximation for the 
effective potential, rather than the random phase approximation (RPA) . In section 5 these 
results were used to obtain an improved estimate for the absorption of radiation by small 
particles, at frequencies small compared to the plasma frequency. 



It is desirable to consider the limitations of some of the approximations which have 
been employed. The RPA prescription itself is an uncontrolled approximation, but it is 
expected to work well at high electron densities (small r s /ao, or equivalently at small values 
of the parameter V defined in (4.21)). 

The image source method simplifies the calculation of the polarisability in two ways. 
Firstly, it treats the surface as if it were flat, and ignores contributions arising from re- 
flections at more distant parts of the surface. These contributions are assumed to have 
relatively smaller amplitude and to combine incoherently, but their effect is very hard to 
quantify (except for slabs or strips with flat parallel faces). A second, and more significant, 
defect of the image source approximation is that it ignores interference effects between the 
direct and reflected paths. These are most significant when the amplitudes for both paths 
are comparable: this happens when either r or r' is close to the surface, and results in the 
polarisability vanishing as either point approaches the surface. The polrisability will differ 
from (2.4) by interference terms which oscillate with a wavenumber comaparble with the 
Fermi wavenumber, kp. 

The semiclasical RPA approximation scheme makes the further assumption that the 
charge density has its classical distribution. In three dimensions the semiclassical RPA 
potential was compared with the image source approximation. The former has a delta 
function at the boundary, whereas the latter has a finite slope at z = 0, corresponding to 
a rapid initial decay with length scale A s = ^e /e 2 is equal to the Thomas-Fermi screening 
length. 

The parameter V which appears in the image source theory for the energy transfer 
on reflection is small when r s /ao is small. In this limit the ratio of the Thomas- Fermi 
screening length to the Fermi wavelength, A s /Af is large. When A s is small compared 
to the Fermi wavelength, the approximations used here are expected to fail, because the 
polarisability approaches zero within Af of the surface, whereas the charge is expected to 
accumulate with a layer of thickness A s . These considerations indicate that the theory is 
more accurate when V is small, which corresponds to the case of good metals with a high 
density of conduction electrons. It is even possible that the results might be asymptotic 
to the results of an exact implementation of the RPA equations in the limit r — > 0. 

In summary, the results presented in this paper represent an approximation scheme 
which leads to analytic expressions for the effective potential (j)(z), the energy transferred 
to an electron rebounding from the surface 8E(9) and the rate of energy absorption by a 
small particle due to an oscilating electric field. The formulae are expected to be a good 
approximation to the full RPA equations when V is small, but it would be desirable to 
benchmark them against a numerical evaluation of the full RPA scheme. 
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Figure captions 

Figure 1. Schematic plot, illustrating the three components of the effective potential. 

Figure 2. The energy 5E transferred to an electron rebounding from the surface, as a 
function of angle of incidence 9. The scale energy is 5E = IqujevvY- This function 
depends upon the material-dependent parameter T oc r s /ao- 



Table caption 



Table 1. Values of the functions ^"(r) and Q(T), defined by (5.25) and (5.30), obtained by 
numerical integration. 



Table 1 
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F(T) 


G(T) 


0.001 


2.679 


0.730 


0.01 


1.594 


0.609 


0.1 


0.807 


0.419 


1.0 


0.813 


0.489 


10.0 


1.091 


0.705 


100.0 


1.149 


0.750 


1000.0 


1.155 


0.755 
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